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C , A modified spin-wave theory is applied to the one-dimensional quantum Heisen- 

berg model with long-range ferromagnetic interactions. 



(TL = —JJ2i<j( r ij) p Si • Sj). Low-temperature properties of this model are investi- 
gated. The susceptibility and the specific heat are calculated; the relation between 
their behaviors and strength of the long-range interactions is obtained. This model in- 
cludes both the Haldane-Shastry model and the nearest-neighbor Heisenberg model; 
the corresponding results in this paper are in agreement with the solutions of both 
the models. It is shown that there exists an ordering transition for 1 < p < 2 where 
the model has longer-range interactions than the HS model. The critical temperature 
is estimated. 
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I. INTRODUCTION 



The quantum Heisenberg chain with long-range interactions decaying as 1/r 2 , which is 
called the Haldane-Shastry model (HS), has been extensively investigated since the exact 
eigenstates and their eigenenergies are obtained independently by Haldane and Shastry. 
The thermodynamics of this model was investigated by Haldane. [[| Unfortunately few 
studies for the general type of l/r p have been made. What we know are rigorous bounds 
for the correlation functions in the disordered phase, which are established by Ito. M 

On the other hand, the model of the classical spins with long-range interactions has been 
studied for about twenty years. In the case of 1-component spins (Ising model), it was proved 
by Dyson || that the model in the region 1 < p < 2 has an ordering transition and that the 
model in the region p > 2 doesn't have it. In the case of 2-component spins (XY model), 
Simanek || showed that with the low-temperature harmonic approximation there exists a 
Kosterlitz-Thouless-like transition to a low-temperature phase with infinite susceptibility. 
In the ferromagnetic case of 3-component spins (Heisenberg model), simulation of the case 
when p = 2 was done by Romano |F| with the Monte Carlo method. The renormalization 
group approach was made by Fisher et al. |8||| 

In comparison with the model with long-range interactions, the Heisenberg model with 
a short-range interaction i.e. a nearest-neighbor interaction (NN Heisenberg model) has a 
longer history. It is well known that for S = 1/2 the exact solution is obtained with the 
Bethe ansatz method. The modified spin-wave approximation was first proposed by one of 
the authors. fH| The conventional spin-wave theory cannot be applied to the case when the 
dimensions are less than three. This modification makes the spin-wave theory valid even 
in one and two dimensions. Results from the modified theory are in good agreement with 
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those from the Bethe ansatz integral equations. []TT1 , |T2 

In this paper we consider the quantum Heisenberg model with long-range ferromagnetic 
interactions decaying as l/r p in one dimension. Its Hamiltonian with a periodic boundary 
condition is written as follows; 

i N N-l JL 

n = -o E £ J l—?^rs m ■ s m+n (i.i) 

z m=\ n=l blIi V AT J 

This model in the limit p — > oo is the NN Heisenberg model and this model for p = 2 
is the HS model. We discuss a modified spin-wave theory of this model and study the 
properties at low temperature. In the next section a formulation of the modified spin- 
wave theory is done. In §3 the terms up to the second order of Bose operators in the 
transformed Hamiltonian are considered. In the region p > 2 the temperature-dependence 
of the susceptibility and of the specific heat are calculated; the critical temperature in the 
region 1 < p < 2 is estimated. In §4 the terms up to the forth order of operators are 
considered. The procedure to obtain physical quantities is shown. The susceptibility and 
the specific heat are analytically obtained for p = 2. In §5 we discuss the results. They are 
compared to the solutions of the HS model and of the NN Heisenberg model. 

II. FORMULATION OF MODIFIED SPIN- WAVE THEORY 

First the Holstein-Primakoff transformation 

$m = + = V^S* f m (S)a m 

S m = — iS"^ = \/2Sa\ n f m (S) 

S m S d m Qj m 



US) = ^l-(-^)aU» = 1 - ^ala m + 0(S- 2 ) (2.1) 



is applied to the Hamiltonian ( |1 . 1|) . Expanded with respect to 1/S, the Hamiltonian is 
rewritten as follows: 

n = E + n 2 + n i + o(s~ 1 ), (2.2) 

-I N N-l TT_ 

ft = -5 E £ ^idfef. 

z m=l n =l blI1 v JV / 

J W iV-1 7T, 

^2 = - ^2 X! • /^n \] P ( a m+n ~ 0, m )( a m+n — a m ), 

/ m=l n=l S1I H AT ) 

j JV AT-1 7L 

= q /7m\ ] P { a m a m+n( a m — a m+n) 

m =l n=l S1I H JV ^ 

Next the site representation is changed to the momentum representation with the Fourier 
transformation a m = (l/VN)J2 k e ikm a k , = (l/^N) J2 k e~ ikm a\. Then H 2 and ft 4 are 
transformed to the following equations: 

N-l TT_ 

k n=l S1I H N I 

+ri(k l + k 2 - k 3 ) - 2r](k 1 - k 3 ) - rj(k 3 - k 2 ) - r)(k 2 - h)}. 

where Re[ x ] stands for the real part of x. 

We consider the expectation value of 7i for the state 

iK}>=nKo^(4no> (2.3) 

k 

which is an eigenstate of 7i 2 . Then we get 

E = E + < U 2 > + < Hi >, (2.4) 
< U 2 > = Y,<n k > S{ V (0) - Re[ V (k)}}, 



<n*> = -^(E < nl > 2 < nk > {Mo) - m - v(-k)} 

k 

+ <n k n k ,>{r)(0)+r)(k-k')-r){k)-r)(k')}), 
k,k'{k^k') 

where n k = a\a k . Magnetization in the z-direction is given by SN — Y,k a \ a k, so the zero- 
magnetization condition is 

SW-5>la fc = 0. (2.5) 

k 

We assume here that n k is the Bose distribution; the entropy and the free energy are 
respectively written as follows: 

(entropy) = ^{(1 + n k ) ln(l + n k ) -n k \nn k }, (2.6) 

k 

F = E-Tx (entropy). (2.7) 

Moreover < n\ > which expresses the expectation value of n\ is 2n\ + h k in terms of n k 
(=< n k >). We want to know h k which minimizes the free energy under the constraint 
condition of zero magnetization, so we introduce the Lagrange multiplier ji and minimize 
the following quantity W: 

W = F-nC£n k -SN). (2.8) 

k 

From dW/dn k = the Bose-Einstein distribution 

hk = e/»fr(*)-M) - 1 (2 ' 9) 

is reproduced where e(k) = dE/dh k and f3 = T~ l . From dW/dfi = we have the self- 
consistent condition 



S = ^X> (2-10) 



which determines the chemical potential /i. 

Using the rotational averaging, we obtain the static susceptibility; 



X=^Y,(4+n k ). (2.11) 

k 



III. QUADRATIC THEORY 



In this section we consider the first two terms up to the quadratic term of operators 
in the Holstein-Primakoff transformed Hamiltonian ( |2.2| ). Then the dispersion relation is 
written as follows: 



JV-l 7T 

N lpj 



<k) = JSj2[-^—ni-cos(kn)}. 



n=r sm (f ) 



(3.1) 



Here we take the thermodynamic limit iV — > oo; the dispersion relation is rewritten to 

1 — cos(kn) 



e{k) = 2JSJ2 



n=l 



(3-2) 



The Bose-Einstein integral function gives us the dominant term of the dispersion for any 
positive and small k for arbitrary p (> 1). 



e{k) 



JS((p-2)k 2 (p>3) 
-JSk 2 Ink (p = 3) 
JSuj(p)k p - 1 (1 < p < 3) 



u(p) 



71 



r(p) cos[tt(p- 2)/2] 



(3-3) 
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In the region p > 2, we can determine the chemical potential \x from the self-consistent 
condition Q2.10| ). We have no Bose condensation which breaks this condition ( |2.10| ). The 



satisfaction of this condition means that the system has no ordering transition. We use the 
determined fi to calculate the susceptibility and the specific heat at low temperature for 
p > 2. The continuum approximation of the state density is valid in the region p > 2; then 



the self-consistent condition (|2.10|) is rewritten to 



~~ Ti Jo efeW+v - 1 ~~ 7T Jo de eP £+v - 1 ^ ' ' 

where v = —/3fjb. 

First we consider the region p > 3. From the dispersion relation and the self-consistent 
condition (|3.4|) : 

1 e~itfe 

O = 



2n^/JS((p-2) J o 

1 r°° x~^dx 



2ttJ(3JSC( P -2) Jo t x+v - 1 



1 



2*pJSC{p-2) 
we obtain v as follows: 



n\)(n\)v^+a\) + ■■■), (3.5) 



v- 1 = AC(p-2)S 3 pj (3.6) 

where the Bose-Einstein integral function is used. From eq. ( [2.1 1| ) the susceptibility in this 
region is calculated as follows: 

, K wg-v^ (37) 

In the same way, we have 



(S p [K(p)} p - l uj(p)(3J)~ (2 < p < 3) 
exp{(3JS 2 n 2 ) (p = 2) 

where n{p) = (p — 1) sin[7r/ {p — 1)]; then we have 

^Mp)]^ {S*k{ P )(3J}^ (2<p<3) 
^exp(/?JSV) (p = 2) 
In this section the free energy per site is given as follows: 

F E 1 T 
" TV = ~N + N^ hk£ ^ - j^Y.ii 1 + ^H 1 + n k ) - n k \nfi k } 

k k 



x = < 



(3.8) 



In the limit N — > oo, we have 



f = eo + Sfi + - — In 1 

7T JO (fc 



(3.9) 



where E /N — > eo ( = —JS((p) ) as iV — > oo. From the dispersion relation the free energy 
is obtained as follows: 



f-e p 
T 



~ < 



-7r- 1 A(p)[5a;(p)]^(^)FT (1< p < 3) 



where A(p) = r(l/(p— l))c(p/(P — 1))/G° — 1)- Then we can calculate the specific heat per 
site; we have 



3C(3/2) 
4 v /27rC(p-2) 



? (P>3) 



^ TF A(p)[o;(p)/2]^(5)^ (K P <3) 
for S 1 = 1/2. 

On the other hand, we cannot determine the chemical potential \x in the region 1 < 
p < 2 as can do in the region p > 2. This is because the zero-magnetization condition is 
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broken by the Bose condensation. Instead we can know the critical temperature T c (=l/j3 c ) 
of spontaneous magnetization by the estimation of the critical temperature of the Bose 
condensation. The critical temperature in the region where p is less than 2 and where p is 
in the vicinity of 2 can be obtained from the constant terms of the Bose-Einstein integral 
function. We should use e{k) = JS{u(p)k p ~ 1 + ((p — 2)k 2 } as the dispersion relation to 
obtain a better estimation of T c near p = 2; then we have 

T^JSuipMp)}?- 1 (3.10) 

7lS(p~ 1) 



We note that e(k) — > e(7r) for any k > 2 as p — > 1+. The critical temperature for p which 
is greater than 1 and which is in the vicinity of 1 can be estimated as follows; 

(2 P - 1)C(P) 



2P~ 1 ln3 

for S = 1/2. 



-J (3.11) 



IV. QUARTIC THEORY 

In this section we consider the first three terms up to the quartic term of operators in 



the Hamiltonian (2.2). Because the dispersion relation is given by e(k) = dE/dn^, we have 



N-l JT_ 1 1 

< k ) = E J [— feH 1 - ™s(kn)}iS --E^' + ^£n fc ,cos(£/n)]. (4.1) 

n=l S1I1 l AT J iV fc' iV k' 

The even function in respect of k is expanded into the following Fourier series: 



f(0) 00 

n k = J -Y + E /H cos(fcm), (4.2) 

Z m=l 

f{m)= { — hkcos(km), (m = 0, 1, 2, ■ ■ •). (4.3) 

•/-7T 7T 

The dispersion in the thermodynamic limit N — > 00 is expressed by 

£ ( k ) = J E ^-t-ffa) ( 4 - 4 ) 

n=l n 

where we use the self-consistent condition of zero magnetization ( |2.10| ) in this limit i.e. 
S=f(0)/2. 

So the problem of calculating physical quantities is reduced to obtaining the dispersion 
and the distribution which satisfy the three eqs. (J2.9|), ([4.3[) and (f4.4|) under the self- 



consistent condition of zero magnetization. For arbitrary p (> 2) we can obtain e(k) and v 
by the following equations: 

= r fl .. f * . (4.6) 



The dispersion and chemical potential which are obtained from the iteration of eqs. fl4.5|) 
and (|4.6j ) give us physical quantities. 

Fortunately we can make an analytical treatment for p = 2; it is shown in the rest of 
this section. Equation (JO) for p = 2 is differentiated twice; we have 

% = Jnu~ JS. (4.7) 

The introduction of a function g(k) = de/dk and the integration of the differentiating 
equation ( |4.7| ) give us the following equation: 
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2 2J 
q = — In 



1 - e~P £ - v 



- 2JSe (4.8) 



1 - e~ v 

where the initial conditions g = e = 0aik = are used. From eq. ( |4.8|) , we use the 
equation: 

de 

dk = " (4.9) 

' 2J ln|^^|^2J<^ 



/3 I l-e-° 

to change an integral variable from to e in the integrations 
The substitution of tt for k in eqs. (|4.4|) and (|4.8[) gives 



e(7r) = ^ + 0(«lni), (4.10) 

e" = 1 - ; - TT (4.11) 

respectively. At low temperature, then we have 

w~exp(-— — ). (4.12) 

From eqs. (|4.9| ) and ( p.llj ) the susceptibility at low temperature of the HS model is 

calculated as x — (P / 6) \J2 / P Jn exp(P J S 2 n 2 / 2) . For S = 1/2 the susceptibility is obtained 
as follows: 



P [2 pjn 2 

x ^6iw^ eM ^ ] - (413) 

The free energy per site in this section is given as follows: 

1 OO T OO f( \ 

/ = iE^{(/(«)} 2 -^E^+^ 

^ n=l '* n=l '* 

1 reM die 

+J_ r >a ^\n{l-e-^- v )de. (4.14) 
pix Jo de 

Then we obtain the following dominant term of the specific heat per site: 

c^l(j) (4.15) 

for p = 2 and 5 = 1/2. 
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V. DISCUSSION 



In this paper we have investigated the one- dimensional quantum Heisenberg model with 
long-range ferromagnetic interactions by the modified spin-wave approximation. This ap- 
proximation makes it possible to treat the cases not only of the special values of p. 

We can apply the modified spin-wave theory to the NN Heisenberg model; its validity 



has already been checked. [JTOJ] The limit p — > oo shifts both the quadratic and quartic theory 
in this paper straightforward to the cases of the NN Heisenberg model. 

The susceptibility and the entropy density of the ferromagnetic HS model are obtained 
by Haldane 0] as follows: 



X = f r /2 dvexp(-2(3Jh) * f\/^exp(^), (5.1) 

2,71 JO 4 V P J7T 8 



s=-f dv{\n[2cosh(f3Jh)] - /3Jhtfmh(3Jh)} (5.2) 
n Jo 

respectively, where h(v) = [v 2 — (tt/2) 2 }/A. From eq. ( |5.2j ) the specific heat per site is 
obtained; we have 

*=!<£ + •- («) 
Both the results of the specific heat in the quadratic theory and that in the quartic theory 
are the same as this specific heat (|5.3|) . Our result of the susceptibility from the quadratic 
theory is not the same as the expression ( |5.1|) but is in agreement with it from the point of 
view that both have exponential divergence at low temperature. So it is reasonable that the 
quadratic theory is qualitatively valid. 

From our results in the quadratic theory we can divide the region of p which determines 
the strength of long-range interactions into two parts. One is the region where temperature- 
dependence except for coefficients of physical quantities is the same as the case of the NN 
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Heisenberg model at low temperature. The other is the region where at low temperature 
it is different from the case of the NN Heisenberg model. We will call the former region 
the effectively short-range region (ESRR) and we will call the latter region the essentially 
long-range region (ELRR). In the region p > 3 the susceptibility has the same power of T 
as that of the NN Heisenberg model in spite that the model in this region has long-range 
interactions. The specific heat also has the same power of T. In the region 2 < p < 3 the 
susceptibility and the specific heat have behaviors of power of T as temperature goes to 
zero. The powers of temperature in the expressions of the susceptibility and of the specific 
heat, however, are changed as p goes from 3 to 2. Especially the power of the susceptibility 
goes to infinity as p — > 2+; the susceptibility for p = 2 has exponential divergence at low 
temperature. Then the region p > 3 is the ESRR and the region p < 3 is the ELRR. The 
authors believe that the model for p = 3 is in the ELRR. This is because there is a possibility 
that the susceptibility has the divergence of 1/T 2 with a factor of correction ln(l/T). The 
conclusion of the quadratic theory is summarized in Fig. 1. 

Until now no approaches to the estimation of the critical temperature in the region 
1 < p < 2 have been known. So we cannot compare our results with others. Our estimation 
of T c , however, is convincing from the following two points of view. One is that T c of the 
expression (|3.10|) vanishes as p — > 2. The other is that T c of the expression ( [3.1 1|) is divergent 
as p — ► 1. This has no contradiction to the fact that it takes an infinite amount of energy 
to make one-magnon state from the vacuum. 

Our result ( }4.13|) of the susceptibility in the quartic theory are in good agreement with 
the expression (|5.1| ); the only difference is the constant factor 2/3. This constant factor 
comes from the point that the modified spin- wave theory is not rotationally invariant. This 
constant factor reminds us of the difference between the susceptibility from the modified 
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spin-wave theory and that from the Schwinger boson mean field theory. ]T3[ Results from 



the Schwinger boson mean field theory will be published elsewhere. There it will be shown 
that almost all the same discussion in this paper is made again. The only difference will be 
the improvement of this constant factor in the susceptibility. 
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is often used in this paper. Analytical property of this function near v = is known; we 
have 



APPENDIX A 



The Bose-Einstein integral function 




(A.l) 



F(a,v) = T(l-a)v a - 1 + 



oo 



((a — n) 
n! 




a—1 a—1 -i 



n! 



F(a,v) 



_{£--!»„} + £ 

L J- r=l ' n^a- 



(-v) n (aeN), 



where C( a ) is Riemann's zeta function. We need the summation 



E 



oo 



1 — cos(kn) 



(A.2) 



n=l 
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to obtain the dispersion relation. We can calculate the dominant term of this summation 
for small k using F(a,v); ( |A.2 ) is expressed by Ci a ) ~ Re[F(a, —ik)]. The following 



integrations are very useful to obtain the susceptibility and the free energy: 



e x+v _ ;p 



T(a)F(a-l,v), (A.3) 



/•CO 

/ x a - 1 \n{l-e~ x ' v )dx = -r(a)F(a + l,v) (A.4) 

JO 



respectively. 
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Figure Captions 



Fig.l The critical temperature in the region 1 < p < 2 and the exponent of the 
susceptibility in the region p > 2. The solid line stands for the exponent. The exponent for 
p = 3 isn't obtained in this paper. The critical temperature is expressed by the broken line. 
Numerical results link the analytical results (|3.10|) and fl3.11|) . 
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